Quadrupole Oscillation of a Single- Vortex Condensate: Evidence for Kelvin Modes 
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We study the two transverse quadrupole modes of a cigar-shaped Bose- Einstein condensate with 
a single centered vortex. We show that the counter-rotating mode is more strongly damped than 
in the absence of a vortex, whereas the co-rotating mode is not affected appreciably by the vortex. 
We interpret this result as a decay of the counter-rotating quadrupole mode into two excitations of 
the vortex line, the so-called Kelvin modes. This is supported by direct observation of the wiggling 
vortex line. 
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When a superfluid described by a macroscopic wave 
function ip(r) is set into rotation, quantized vortex lines 
appear, along which the density \ip\ 2 is zero [jj], |[ |J. 
Since ip is single valued, its phase variation on a closed 
contour around a vortex line is 2mr (n is integer). Such 
vortices have been observed in many systems, e. g. super- 
conductors [Q, superfluid liquid helium j^], and gaseous 
Bose-Einstein condensates q, 0, [|, || . 

A vortex line is a dynamic object. As for a classical 
string, transverse vibration modes of the line can be ex- 
cited. For a classical vortex, the modes have been calcu- 
lated by Kelvin in 1880 0. His result can be transposed 
to a single quasi-linear quantum vortex with unit charge 
(n = 1) in a homogeneous superfluid. This gives the rela- 
tion between the energy Hu>k of a quantum of the Kelvin 
modes (kelvon) and its wave vector k [[| [uj [l2| : 



^ln(l/feO (*£<!)• (1) 

Here m is the mass of a particle of the fluid and £ = 
{Sirpa) -1 / 2 the healing length (a is the scattering length 
characterizing the binary interactions in the fluid and p 
the fluid density). Due to the Kelvin-Helmoltz theorem 
[ tuj , the Kelvin modes rotate always in the sense oppo- 
site to the vortex velocity field || |l0) . Consequently the 
angular momentum of a kelvon associated with a n — 1 
vortex is —h [0. 

The Kelvin modes play an important role in the dy- 
namics of superfluids || 15 or neutron stars ]l3| . In 
superfluid liquid helium 4 He, they have been observed 
by trapping ions in the vortex core, and exciting them 
by a circularly polarized electromagnetic field Q . 

In this Letter, we present experimental evidence for 
Kelvin modes excited through non-linear (Beliaev) de- 
cay of the quadrupole mode m = —2. An elementary 
excitation of the quadrupole mode m = —2 decays into a 
pair of kelvons with wave vectors k and —k. We measure 
the corresponding increase of the decay rate with respect 
to the decay of the m = +2 mode, for which angular 
momentum conservation forbids such a mechanism. We 
also show pictures of the wiggling vortex line obtained 
after the excitation of the m = —2 mode. The spatial 



period of the oscillations of the line is in good agreement 
with the one deduced from the generalization of Eq. (|l|) 
to a trapped condensate [J17| [l8| and the energy conser- 
vation = 2u>k (uj-2 is the frequency of the m = —2 
quadrupole mode). 

We use a cigar-shaped 87 Rb condensate held in an axi- 
symmetric Ioffe-Pritchard magnetic trap. The atoms are 
spin-polarized in the F = mp = 2 state. The trapping 
frequencies are lu z /2tt =11.8 Hz (longitudinal axis) and 
oj±/2tt = 98.5 Hz (transverse plane). The temperature 
is ~ 70 nK, corresponding to a 70% condensed fraction. 
The number of atoms in the condensate is 1.3 x 10 5 and 
the chemical potential is 40 nK. 

In order to produce in a reliable way a single vortex 
centered on the z axis jl9| , we follow the procedure out- 
lined in Ref. j2(j]. We use an off-resonant laser beam, 
whose motion is controlled by acousto-optic modulators, 
to superimpose a rotating dipole potential onto the mag- 
netic potential. This dipole potential is asymmetric in 
the x—y plane and reads muj 2 L e(X 2 —Y 2 )/2 with e = 0.07. 
The axes X, Y are deduced from the fixed axes x,y by a 
rotation of an angle fit. The dipole potential is switched 
on after the condensate formation for a period of 0.3 s. 
We choose Q close to the rotating quadrupole resonance 
w_l/v2j thus creating a lattice with 5 to 7 vortices. Dur- 
ing the next 2 seconds, the laser stirrer is blocked and 
the vortex lattice decays due to a slight (~ 1%) static 
anisotropy of the magnetic potential in the x — y plane. 
The condensate is then left with a single centered vortex. 
The lifetime of this last vortex is ~ 7 seconds, which is 
much longer than the rest of the experimental sequence. 
Experiments without vortices performed for comparison 
follow the same procedure, except that rotation is kept at 
a lower frequency to prevent any nucleation of vortices. 

The dipole potential is then used again on the sin- 
gle vortex condensate to selectively excite the modes 
m = +2, m = —2 or their superposition. During this 
excitation, the angular frequency and deformation pa- 
rameter of the dipole trap are set to ft' and e'. The du- 
ration of the excitation is denoted by r. The state of the 
condensate is finally probed after a 25 ms time of flight 
by simultaneous absorption imaging along the two direc- 
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FIG. 1: Angle 9 and ellipticity ( as a function of t in absence 
(a), and in presence of a single centered vortex (b). In (b) the 
precession of the main axes observed for small t increases at 
longer times, and the jumps in 9 become rounded indicating 
that the mode m = —2 decays faster than the mode m = +2. 



tions y and z. Transverse images (imaging beam along 
y) give indications on the behavior of the vortex line. 
Quantitative information is obtained from the longitu- 
dinal images (imaging beam along z). These are fitted 
assuming an elliptic Thomas-Fermi profile for the spatial 
density of the condensate in the x — y plane. We mea- 
sure the size R\ and the polar angle 9 of the long axis, 
and the size R s of the short axis. C = R\/Rs denotes the 
ellipticity in the x — y plane. 

The first series of experiments aims at observing the 
free evolution of the two quadrupole surface modes to = 
±2, in the presence and absence of a vortex. A per- 
cussional excitation is performed using the laser stirrer 
with fixed axes (fi' = 0) for a short duration r = 0.5 ms 
<C ujJ 1 (e' ~ 1). This excites a superposition of to = +2 
and m = —2 modes with equal amplitudes. We then let 
the cloud evolve freely in the magnetic trap for a variable 
time t and we perform the time of flight analysis. The 
quantities 9 and £ are plotted as a function of t in Fig. [l]. 

In the case where the condensate is vortex-free 
(Fig. |I|a), the angle 9 jumps periodically between 
and it/ 2, indicating that the amplitudes of the two 
quadrupole modes m = +2 and to = —2 stay equal. The 
situation is dramatically different in the presence of a vor- 



tex (Fig. |l|b). In this case two successive regimes occur: 
for short times, one still observes the quadrupole oscilla- 
tion, now with precessing axes. The precession is due to 
the lift of degeneracy between the frequencies uj±2 of the 
two modes m = ±2 with uj +2 > w_ 2 © ||, |§ |§. This 
well known effect is often used to measure the angular 
momentum of the condensate |2^, [27], ^8|. For t > 5 ms, 
the precession rate increases and the jumps of tt/2 in 9 
become more and more rounded. This behavior indicates 
that the mode m = —2 decays faster than m = +2, which 
will eventually lead to an atom cloud with non-oscillating 
ellipticity and constant rotation at uj + 2- 

For a quantitative analysis, we have fitted 9(t) and £(i) 
assuming that the amplitudes a m of the modes m = ±2 
vary as 



+ («w m + r,„) a m = o 



(2) 



with a + 2(0) = a_2(0). The decay of the quadrupole 
modes is modelled using phenomenological linear damp- 
ing rates T m . A more refined treatment should take into 
account the non-linear character of the Beliaev-type pro- 
cess under investigation. 

In the case of no vortex (Fig. |l|a), where the sym- 
metry between the modes m = ±2 is preserved, we 
deduce from the fit J^^tt = 142.0 ± 0.5 Hz, which 
is close to the prediction in the Thomas-Fermi regime 
V2w_l/2tt = 139.3 Hz Q. For the decay rate we 
find = 21.3 ± 1.3s -1 . In the presence of a vor- 
tex (Fig. |l]b), the fit gives: lu +2 /2tt = 159.5 ± 1.0Hz, 
w_ 2 /2tt = 116.8 ± 1.0 Hz, r +2 = 19.1 ± 2.0s- 1 and 
T_2 = 35.7±4.0s -1 . We note that the two measured fre- 
quencies satisfy the sum rule uj\ 2 + ^-2 = 4wj p5| with 
a good accuracy. The fit confirms the difference in the 
two decay rates with the addition that the mode m = — 2 
decays faster in the presence of a vortex (T_ 2 > F±l), 
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whereas T 

In a second series of experiments, we perform a 
spectroscopy-like study of the two modes m = +2 and 
to = —2 separately. For the excitation, we set e' = 0.025 
and apply the laser stirrer for t = 40 ms. It is rotating 
with fl' either with the same (excitation of to = +2) or 
the opposite sense (m = —2) with respect to the vortex. 
We perform the usual time of flight and imaging, and 
measure the ellipticity £ as a function of O' . We observe 
a clear resonance for each mode (Fig. ^). The two reso- 
nances occur at different central frequencies according to 
the already mentioned lift of degeneracy. The key feature 
of the present study is that the width of the m = —2 reso- 
nance is significantly larger than that of to = +2. Hence, 
as seen before, the vortex causes larger damping for the 
counter-rotating surface wave than for the co-rotating. 

To analyze our data, we assume that the amplitude of 
the relevant mode varies according to Eq. (||) with a driv- 
ing term A e~ 2in * on the right hand side. This assump- 
tion leads to a good fit to the data (Fig. @) and yields 
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2: Resonance of the ellipticity £ for the m = +2 and 
m — — 2 modes. The lines are fits deduced from the solution 
of Eq. ^ with a sinusoidal drive. The resonance of the mode 
m — —2 is broader than that of the m = +2, as a consequence 
of a larger damping rate. 



cj +2 /27r = 161.0±1.0Hz, cj_ 2 /27t = 119.8±1.0Hz. These 
frequencies are in excellent agreement with those ob- 
tained by the percussional excitation. The ~ 2 Hz shift is 
due to the increase of the trapping frequency by the laser 
stirrer. The fitted damping rates are r + 2 = 24 ± 5s _1 
and r_2 = 57±10s _1 . These are slightly larger than 
those obtained from the percussional experiment, possi- 
bly due to an additional heating caused by the stirrer. 
However, the main feature T_2 ~ 2 T+2 remains valid. 

We now discuss the possible physical origin of this in- 
creased damping of the to = — 2 mode. The starting 
point of our analysis is the behavior of the uncondensed 
part of the gas. Is it still rotating when we perform the 
percussional or spectroscopic study? If it rotates, the 
symmetry between the two senses of rotation is broken 
and the co-rotating mode experiences less friction than 
the counter-rotating mode. In [p0|, it has been shown 



that in this case T_2 > ^±1 an d T + 2 < T 



,(o) 

±2- 



Moreover 



one gets r_2 + T + 2 — ^±1- This latter relation does not 
correspond to our observations. Although we do observe 
an increase of T_2 when a vortex has been nucleated, we 
do not observe a corresponding reduction of T + 2 ■ From 
this we infer that the uncondensed fraction does not ro- 
tate significantly in our experiment. Furthermore, no 
such rotation is expected for our experimental conditions. 
The ~ 1% static trap anisotropy, which barely affects the 
behavior of the condensate, rapidly damps the rotation 
of a non-condensed gas |31, p2| . The thermal cloud thus 



experiences two competitive forces in the bare magnetic 
trap: a rotational drive from the velocity flow of the con- 
densate, and a friction from the static trap anisotropy. 
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FIG. 3: Transverse images of the condensate after a 33 ms 
resonant excitation of the m = —2 mode (a) and the m = +2 
mode (b). The picture on the left shows a periodic structure 
superimposed on the vortex line. The horizontal density pro- 
files (c,d) have been obtained by averaging over three vertical 
pixels around the center of the vortex line. 



Unlike in the case of a complete vortex array Q, the 
rotational drive due to a condensate with a single vor- 
tex is small compared to the friction term. Therefore the 
rotation of the thermal cloud is expected to be negligible. 

Assuming now that the uncondensed gas is at rest, the 
other possible reason for a faster damping of the to = — 2 
mode is the existence of a decay channel open for to = — 2 
and closed for to = +2. The natural candidate for this 
mechanism is the Beliaev conversion of the to — —2 mode 
into two Kelvin modes k and —k, followed by the kelvon 
decay into another kelvon with a slightly smaller momen- 
tum and a phonon. Since the Kelvin excitations have 
negative angular momentum with respect to the vortex 
charge, this mechanism is specific to the to — —2 mode. 
A recent theoretical analysis of this decay process, taking 
into account both Beliaev and Landau damping, has led 
to damping rates in good agreement with those that we 
measure experimentally pjj . 

To support the interpretation of decay into Kelvin 
modes, we present time-of-flight images of the vortex line. 
The images are taken immediately after the excitation of 
the to = —2 or m = +2 mode by the spectroscopy-like 
method. For each mode we choose f2' at the center of 
the resonance and apply the probe for r = 33 ms. Fig. ^ 
shows two condensates in the transverse view. The left 
picture (a) was taken counter-rotating (to = —2), the 
right picture (b) co-rotating (m = +2). Below are shown 
the horizontal density profiles taken at the center of the 
vortex line (c,d). 

The vortex line is visible in both images with good 
contrast. It has the form of an unfolded "N" pcfl . Af- 
ter excitation of the m = — 2 mode, the atom density 
presents a periodic structure apparent as vertical stripes 
in the vicinity of the vortex line (Fig. ^a,c). A Fourier 
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analysis of the density profile reveals a peak centered 
at ko — 0.70 //m with a full width at half maximum 
Afc — 0.15 /Ltm -1 . The area of the peak, which can 
be observed reproducibly for similar experimental con- 
ditions, is 6 standard deviations above noise. On the 
contrary, such a periodic structure never appears in the 
images corresponding to m = +2 (see e.g. Fig. ||b,d). 

This observation supports the proposed decay mecha- 
nism of the m = — 2 mode into a pair of Kelvin modes. 
At any given time, the excited vortex line is expected to 
have an almost sinusoidal shape with a spatial period of 
2-7r/fc along the z axis. It is located in a plane contain- 
ing the z axis and rotating around this axis in the sense 
opposite to the vortex. At the center of the condensate, 
the healing length is £ ~ 0.3 /im. Using Eq. ([[]) and 
energy conservation, we find k — 0.8 /im -1 . Assuming 
that during the time of flight, the expansion factor A 2 
along the z axis is the same for the shape of the vortex 
line as it is for the longitudinal size of the condensate 
(A z ^1.3 j}5|), we recover the measured wave vector fc 
within Afco. This gives ~ 9 oscillation periods over the 
L ~ 70 /im length of the condensate, which is consistent 
with the observation in Fig. ||a. With the quantization 
condition k = Afir/L for the vortex line, we get Af ~ 18. 
The splitting between the adjacent Kelvin modes is then 
~ 20 s _1 , so that we have one or two Kelvin modes within 
the width T±l of the quadrupole resonance. We note that 
a detailed treatment of the relation between the observed 
density modulation and the oscillations of the vortex line 
remains to be done. The principle of such a treatment 
may follow that of Ref. |36| , where phase fluctuations due 
to a short coherence length in a quasicondensate appear 
as density fluctuations after a time of flight. 

In conclusion, we have observed a strong difference be- 
tween the damping rates of the co-rotating (m = +2) and 
counter-rotating (m = —2) quadrupole modes of a single 
vortex condensate. We explain this difference as due to 
the decay of the m — —2 mode into Kelvin excitations of 
the vortex line. This is confirmed by images of the line, 
which show wiggling only for counter-rotating excitation. 
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